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The electron-electron scattering rate of single particle ex- 
citations in atoms is estimated and compared with the cor- 
responding rate in quantum dots. It is found that in alkali 
atoms single particle excitations do not acquire a width due 
to electron-electron interaction, while in complex atoms they 
may. This width is typically smaller than the single particle 
level spacing, and hence does not affect the number of dis- 
crete single particle excitations resolved below the ionization 
threshold. This situation is contrasted with that of quantum 
dots where electron-electron interaction severely limits the 
number of resolved excitations. Unlike the case of quantum 
dots, the scattering rate in atoms is found to decrease with 
increasing excitation energy. The different effect of electron- 
electron interaction on the spectrum of quantum dots and 
atoms is traced to the different confining potentials in the 
two systems. 



I. INTRODUCTION 

This work investigates theoretically the lifetime of sin- 
gle particle excitations due to electron-electron (e-e) in- 
teraction in atoms. In particular we are interested in 
the manifestation of this lifetime in the measured spec- 
trum. Although a large atom and a quantum dot (QD) 
may contain a similar number of electrons, their spec- 
tral characteristics are very different. In QDs the level 
broadening by e-e interaction severely limits the number 
of discrete single particle levels. For example, Sivan et 
al. Q were able to resolve only ~ 10 discrete levels in the 
excitation spectrum of a diffusive QD containing about 
4000 electrons. In atoms on the other hand, hundreds 
of spectral lines have been measured and tabulated, with 
no appreciable effect of e-e interaction on level widths. 
Hence, the effect of e-e interaction seems to be very dif- 
ferent in these two systems, and this work attempts to 
understand the roots of this difference. We do not intend 
to obtain exact results for a specific atom, but rather to 
develop a general understanding of how atomic spectra 
are influenced by e-e interaction in comparison with QDs. 

Atoms and QDs both contain a comparable number of 
interacting electrons in a confined volume. However, they 



differ in size, symmetry and confining potential. The lat- 
ter are typically larger, possess no particular symmetry, 
and are confined by a sharper potential. We find that 
the different spectral characteristics of these two systems 
are primarily due to the different confining potentials. 

We focus our attention on excited atomic states in 
which one electron is excited to a weakly bound state, 
while the other electrons remain in their ground state 
configuration. The excited electron then occupies a 
hydrogen-like orbital, whose spatial extent is much larger 
than that of the remaining electrons. The excited elec- 
tron is subject to the Coulomb potential of the rest of 
the atom (which we refer to as "ion"). This interaction 
is divided into a static potential (averaged over ionic de- 
grees of freedom) and a residual interaction. To calculate 
the lifetime of an excited atomic state, we consider the 
residual interaction acting on the excited electron in the 
confining potential of the ion. 

The ionic spectrum of alkali atoms is characterized by 
a large gap (on the order of the ionization energy) at the 
ground state. As a result, single particle excitations are, 
to good approximation, eigenstates of the atom, and no 
level broadening due to e-e interaction is expected. Com- 
plex atoms with a number of valence electrons are differ- 
ent. Typically, there are several open shells with many 
single particle states of similar binding energies. The cor- 
responding number of many-body states, N to t, is hence 
exponentially large. We find that the interaction matrix 
elements are cx l/y/Ntot, while the relevant density of 
states is cx Ntot- Consequently, for large Ntot the matrix 
elements of the residual interaction are larger than the 
level spacing of the spectrum to which the excited state 
is coupled. The use of Fermi's golden rule is hence jus- 
tified. The resulting level width is independent of Ntot, 
and typically smaller than the single-particle level spac- 
ing @. 

The different effect of e-e interaction in complex atoms 
and QDs can be qualitatively understood in the follow- 
ing way. The scattering rate is proportional to the resid- 
ual interaction matrix element squared times the den- 
sity of final states (relaxed electron plus excited ion) . In 
QDs, due to the sharp confining potential, the matrix 
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elements are on the average only weakly dependent on 
energy. The density of hnal states is proportional to the 
excitation energy squared. The resulting scattering rates 
grow with energy and eventually exceed the fairly con- 
stant single particle level spacing. In complex atoms, 
due to the shallow confining potential, the radius of the 
excited electron's orbit grows rapidly with energy lead- 
ing to suppressed matrix elements and diminishing e-e 
scattering rates. The decrease of e-e scattering rates as 
a function of excitation energy is in sharp contrast to 
the corresponding trend in QDs and Fermi liquids. The 
hydrogen-like single particle level spacing is also reduced 
with energy, but slower than the scattering rates. Con- 
sequently, in complex atoms the discrete nature of the 
spectrum is preserved. 

The outline of this paper is as follows. In section (||) 
we define e-e lifetime of single particle excitations and 
discuss its meaning in a finite system. We consider in 
particular the manifestation of e-e lifetime in a spectro- 
scopic measurement. In section ([II) we estimate the e-e 
lifetime of single particle excitations in atoms. We distin- 
guish between alkali atoms for which we find no broaden- 
ing due to e-e interaction, and complex atoms for which 
such broadening may occur. Finally, in section (IV) we 
compare our results for e-e lifetime in atoms to previous 
results concerning e-e lifetime in QDs. 



II. ELECTRON-ELECTRON LIFETIME AND ITS 
MANIFESTATION IN THE MEASURED 
SPECTRUM 

Consider an isolated, TV-electron system in its ground 
state, \g.s.), and an excited state, 



l*(* = 0)> 



g.s. 



(1) 



where cj and Cj are single particle creation and annihila- 
tion operators respectively. In the absence of interaction, 
\ty(t — 0)) is an eigenstate, and therefore characterized 
by an infinite lifetime. In the presence of interaction, it 
is no longer an eigenstate and decays with time. The e-e 
lifetime, denoted by r ee , is defined as the decay time of 
the initial state, 



\(*(t)\*(t = 0))f~e- 



(2) 



This analysis assumes zero temperature. Generalization 
to finite temperature is straightforward. 

The meaning of e-e lifetime in a system characterized 
by a discrete energy spectrum, such as an atom below 
the ionization threshold or a QD, should be clarified. Let 
{|V> a )} be a basis of exact eigenstates of the system with 
energies {E a }. The initial state can be expressed as a 
superposition of these eigenstates, 



\*(t = 0)) = 4 Cj \g.s.)=J2*M (3) 

a 

5> Q | 2 = i. 



The time evolution of the initial state is given by, 



!*(*)} = E A « e 



iE a t/h 



\4>a), 



(4) 



and the probability to remain in the initial state at time 
t is, 



|(*(t = Q)|*(t))| 



a e iE a t/n 



(5) 



Strictly speaking, as opposed to the case of a continu- 
ous spectrum, here there is no decay of the probability 
to remain in the initial state. This probability oscil- 
lates, and for long enough times it can get arbitrarily 
close to 1 (it can return exactly to 1 if the energy spec- 
trum is commensurate). Nevertheless, as the number of 
dominant terms in (|j) increases, it takes longer for the 
initial state to reconstruct. The typical time scale for 
this reconstruction is the Rabi-time, TR a u — h/A, where 
A is the average many-body level spacing. On shorter 
time scales the probability effectively decays exponen- 
tially with time. Consequently, we define the e-e lifetime 
in a finite system as in (0), keeping in mind the restriction 
on the time scales t, r ee <C TrcM- 

The e-e lifetime of a single particle excitation can be 
estimated using Fermi's golden rule. The full interacting 
Hamiltonian H is divided into an unperturbed Hamilto- 
nian Ho, which may include part of the e-e interaction 
(for example, the Hartree or Hartree-Fock part), and a 
perturbation Hamiltonian Hi nt which includes the rest of 
the e-e interaction. Let {|/ )} be the iV-particle eigen- 
states of Ho with eigenenergies {E®}. The golden rule 
result for the e-e lifetime of a single particle excitation 
|<°) is, 



(0) 



The sum in (||) extends over all N-particle eigenstates of 
Hq other than the initial one. 

Following 0| , we arrange the eigenstates of the unper- 
turbed Hamiltonian Hq in the form of a hierarchal tree 
emanating from the ground state of Hq. The first gener- 
ation includes all eigenstates of Ho which are coupled by 
the interaction to the ground state, in particular single 
particle excitations (as in ([!])) are included. Each gener- 
ation includes all eigenstates of Hq not previously incor- 
porated in the tree which are connected by non- vanishing 
interaction matrix elements to the previous generation. 
The problem of e-e lifetime of single particle excitations 
can be mapped on the problem of Anderson localization 



2 



on the hierarchal tree, described by a Hamiltonian Ha- 
Each eigenstate |/°) of Hq is a site in the hierarchal tree 
with an on-site energy E°. The hopping amplitude be- 
tween two sites is given by the corresponding residual 
interaction matrix element. The single particle spectrum 
of Ha on the hierarchal tree is equivalent to the many- 
body spectrum of Ho. An excitation of the type (Q) cor- 
responds to a single site on the lattice. Its time evolu- 
tion is determined by the overlap of this site with exact 
eigenstates of Ha- An overlap of this site with extended 
states leads to exponential decay, while overlap limited 
to localized states leads to beating of few frequencies. 

We now turn to discuss how e-e lifetime of single par- 
ticle excitations is manifested in the measured many- 
electron spectrum of a quantum dot or an atom. A fi- 
nite, isolated, Fermi system has discrete energy levels of 
zero width, which correspond to exact many-body eigen- 
states of the system. A spectroscopic measurement in- 
volves coupling the system to some external measuring 
device. Typically, the measurement operators are single 
electron ones. As a result, only many-body eigenstates 
with finite overlap with single particle excitations can be 
detected. 

Consider a specific example, an optical absorption ex- 
periment on a many-electron system in its ground state 
\g.s.). The measurement operator in this case is given 
by, 

= ^Aycjc,- + fc.c., ( 7 ) 

where Ay are coefficients (including matrix elements of 
photon operators) determining the strength of the cou- 
pling to various single particle excitations. The absorp- 
tion of the system is characterized by the spectral func- 
tion, 

«M = E K/l Eij Aycta + h.c.\g.s.)\ 2 (8) 
/ 

x 5{E gs +fiio- E f ), 

where |/) and Ef are exact final many-body eigenstates 
and eigenenergies of the system, and E gs is the ground 
state energy of the system. In the absence of e-e interac- 
tion, a(uj) exhibits a series of (5-peaks corresponding to 
single particle excitations of the system. When e-e inter- 
action is included, each single particle excitation becomes 
a superposition of several exact many-body eigenstates. 
The spectral function then displays many more absorp- 
tion peaks corresponding to many-body eigenstates that 
overlap with single particle excitations. A many-body 
eigenstate can overlap several single particle excitations, 
which then generate interfering contributions to its ab- 
sorption peak. 

A significant simplification of the spectral function and 
its relation to the concept of e-e scattering rate is ob- 



tained by assuming that each many-body state |/) over- 
laps with at most one single particle excitation, an ap- 
proximation equivalent to the reduction of the hierarchal 
tree into a Cayley tree Q. This approximation is valid 
when the mira-generation matrix elements are negligible, 
or when the single particle level spacing is much larger 
than the resulting e-e broadening. Within this approx- 
imation, the intensity of peaks having significant over- 
lap with a particular single-particle excitation c\ Q Cj \g.s.) 

is approximated by W j \ 2 \{f\c\ Q Cj a \g.s.)\ 2 . The relative 
intensity of all 5-peaks |/) associated with this single par- 
ticle excitation is proportional to | (f\c] g Cj \g.s.) | 2 . The 
factor |Ai 0j0 | 2 is common to all of them. 

Consider the time evolution of a system initially in a 
state, \tp(t = 0)) = 4 c jo \g.s.), 

\m)=J2 eiEft/n (M ^o\9-s.)\f)- (9) 

The e-e scattering rate is determined by the energy 
spread of exact many-body eigenstates, |/), that partici- 
pate in the sum (||), provided there are many such states. 
These many-body eigenstates are exactly those that gen- 
erate the absorption peaks associated with c\ o Cj \g.s.) in 
the spectral function. Thus, e-e scattering rate of a sin- 
gle particle excitation is manifested in the energy width 
of the cluster of absorption peaks associated with it in 
the spectral function. 

In a real experiment the resolution of the measuring 
device is finite. The measured many-body spectrum is 
smeared, so that each <5-peak in the spectral function ap- 
pears as a broadened peak. This broadening may have 
a typical scale larger than the many-body level spacing. 
In this case, the many-body eigenstates can no longer be 
resolved in the spectrum. Then, a cluster of 5-peaks as- 
sociated with a given single particle excitation appears as 
a broad single particle resonance whose width is equal to 
the inverse e-e lifetime of the excitation. Single particle 
resonances can be resolved as long as their broadening 
is smaller than the single particle level spacing. When 
the broadening of single-particle resonances exceeds the 
single-particle level spacing, the measured spectrum be- 
comes essentially continuous. 

III. E-E LIFETIME IN ATOMS 

In this section we investigate the e-e lifetime of single 
particle excitations in atoms or ions below the ionization 
threshold || • We consider an atom (or ion) composed of 
a fixed nucleus with charge Z and N electrons. The full 
non-relativistic Hamiltonian of such a system is, 
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where pi and r, are the momentum and position opera- 
tors of the electrons, respectively, m is the electron's mass 
and e is the electron's charge. Since an exact solution to 
the full Hamiltonian is generally not known, one resorts 
to approximation methods. A standard approximation 
is the replacement of the full Hamiltonian by a single 
particle Hamiltonian that includes an effective potential 
induced by the average electron density (e.g. the Hartree 
approximation). The resulting single particle spectrum 
is degenerate due to spherical symmetry, leading to the 
well known shell structure of atoms. 

The degeneracy in the single particle spectrum leads 
to a much larger degeneracy in the many-body spectrum. 
For example, the ground state configuration of an Eu- 
ropium atom includes 7 electrons in an open f-shell. The 
14-fold single particle degeneracy of an f-shell, yields a 
14!/(7!) 2 = 3432-fold degeneracy in the non-interacting 
many-body ground state. Within such degenerate sub- 
spaces the interaction matrix elements are in many cases 
comparable or larger than the many-body level spacing. 
As a result, e-e interaction has significant effect on the 
many-body eigenstates and eigenenergies. This is demon- 
strated in the numerical calculations of the eigenstates 
and eigenenergies of a Cerium atom by Flambaum et al. 

Thus, in a perturbative calculation of e-e lifetime 
most of the interaction has to be included in the unper- 
turbed Hamiltonian, Ho, which is no longer a simple sum 
of single particle Hamiltonians. 

We consider the e-e lifetime of high energy single parti- 
cle excitations (but still below the ionization threshold). 
As a result the excited electron spends most of the time 
in regions where the density of the other electrons is ex- 
ponentially small. Two major simplifications can hence 
be made, 

1. The exchange integral between the excited electron 
and the rest of the electrons is exponentially small 
and can be neglected. The excited electron can 
therefore be considered distinguishable from other 
electrons. 

2. The excited electron is subject to a potential which 
is roughly that of a hydrogen atom. Its wavefunc- 
tion may then be approximated by the correspond- 
ing hydrogenic one, except for s-shell electrons that 
penetrate into the ion. 

The atom is divided in our treatment into an excited 
(distinguishable) electron and an ion containing all other 
electrons and the nucleus. The Hilbert space is spanned 
by direct products of the ionic states and the excited 
electron states. 

The unperturbed Hamiltonian is, 

->2 



Hq — Hi OT: 



P 

2m 



V e (r). 



(11) 



Hi on is the full Hamiltonian of the ion, and h e is the ef- 
fective single particle Hamiltonian of the excited electron 



which includes an effective potential V e induced by the 
nucleus and the spherically averaged ionic ground state 
electron density. In regions exterior to the ionic elec- 
tron density, the effective potential experienced by the 
excited electron is simply V e (r) ~ — e 2 /r. The perturba- 
tion Hamiltonian is, 



Hint — ^ 



dr 



1 Sp(r') 



5p(r') = p{r') — p(r') 
1 



(12) 
(13) 



p(r>) - — I dn(g.s\p(r')\g.s.), 



where p(r') is the density operator of the N— 1 electrons 
of the ion, and Sp(r') is its fluctuating part. \g.s) refers 
to the ionic ground state, and r is the position operator 
of the excited electron. 

The golden rule is now employed to the calculation of 
the lifetime of an excited electron due to the residual 
coupling Hint with the rest of the electrons. We first 
consider the relevant density of states and the interaction 
matrix elements. 



A. The Density of States of Ho 

The density of states of Ho, G{E), is a convolution of 
the density of states of the excited electron, g e (e), and 
the density of states of the ion, gi on (t), 



G{E) 



de g e {e)gion(E ~ e), 



(14) 



where all energies are measured relative to the ground 
state. In order to study the spectrum of Hq, we first ex- 
amine the single particle spectrum of the atom within the 
Thomas-Fermi model, which gives a fairly good picture 
of a large enough atom. We later consider the charac- 
teristics of the many-body spectrum of the ion and the 
effect of e-e interaction. 



1. The Single Particle Spectrum within the Thomas-Fermi 
Model 



The Thomas-Fermi model (TF) (see e.g. [J7j) is a self- 
consistent single particle model used to obtain the ap- 
proximate ground state electron density of an atom from 
which an effective single particle potential is calculated. 
Here, we use it to estimate the density of states for exci- 
tations above the atomic ground state. The TF effective 
potential in a neutral atom is the solution of the TF 
equation, 



d 2 V(r) 
dr 2 



8V2 1 

' 3^ 3/ 2 



{p~V{r)f 2 , 



(15) 
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with, 



V(r -> 0) 



Ze 2 



; V(r oo) ^ 0. 



(16) 



Here, jj, is the chemical potential and ao = h 2 /me 2 is the 
Bohr radius. We define, as customary, 



x = JL Z ^b . 6 = 2 (A 
a \3tt 



2/3 



V(r) - fi 



Ze 2 /r 

With these definitions, the TF equation reduces to. 

* l/2 S =x3/2(x) 

X (0) = l ; X (x^oc)=0. 



(17) 
(18) 

(19) 



This non- linear differential equation is solved numerically 
0, and we denote its solution by Xo( x )- Since the equa- 
tion and its boundary conditions are independent of Z, 
so must be its solution xo(x). The TF potential is given 
by, 



V TF (r) 
Asymptotically, 

V TF (r) ~ \ 



Z£_ 

r 



ao 



.^L r«Z- 1 /3 ao 
-±^4# r»Z-V3 ao . 



(20) 



(21) 



The length scale Z~ 1 > 3 oq characterizes screening in a TF 
atom. 

An obvious flaw in the TF potential is that it includes 
the electrostatic self-interaction of the electron. This has 
a noticeable effect at large distances where the electron 
is essentially outside the ionic charge distribution, and 
therefore experiences a potential V(r) ~ —e 2 /r, rather 
than ~ 1/r 4 as obtained from the TF model. The TF 
potential can be modified to correct this erroneous result 
in various ways [pi leading to, 



V TF {r) 



r < Z-Woq 
Z~ 1/3 a < r < a 
r » ao- 



(22) 



The modification to the potential is important for the 
excited electron we consider, since it spends most of its 
time at distances larger than ao from the nucleus. 

The ionization energy of the atom is the difference be- 
tween the chemical potential of a neutral atom (/i = 0), 
and that of an ion (N = Z — 1). To lowest order in 1/Z, 
the ionization energy is constant and given by B, 



E. 



,.,„ ~ 0.109—. 

ao 



(23) 



The ionization energy is approximately the same for all 
atoms. This result does not account for the irregular 
variation of Ei on as a function of Z arising from the shell 
structure. 

The single particle density of states of a TF atom is 
obtained by numerically solving the Schroedinger equa- 
tion with Vtf- Alternatively, the single particle density 
of states can be calculated in the semiclassical approxi- 
mation employed in the derivation of the TF equation. 
Since we are interested only in the approximated density 
of states |l(J, the semiclassical approach is sufficient. 

The density of states of a free 3D Fermi gas enclosed 
in a volume V is given by g(e) — (2m) 3 / 2 / (2ir 2 h 3 )y/e V . 
To obtain the density of states in a TF atom with a 
space dependent potential we use the expression for the 
free Fermi gas locally, replacing e by e — Vtf i r ) a t each 
point. Integrating over space we obtain the single particle 
density of states in the semiclassical approximation, 

(977,13/2 r R(e) 

g TF (e) = 2 [ —l^ r 2 ^/(e-V T F(r))dr. (24) 

nn Jo 

where R(e) is the classically accessible radius given by 
V TF (R(e)) - e. 

Introducing explicit Z-dependence and using the defi- 
nition of x in (|l^) we obtain, 



g T F(t,Z) 

(2m) 3 / 2 f R ^ 



(25) 



- 2 v -' / 3 / r 2 y/(e-V T F(r,Z))dr 
irn Jq 

1 / Q \ 5 /3 
1 / 37T \ „_ 1/3 0o 



* It z 



XUa /Z i/3 e 2 ) I TT 

2 / ea , X0{X) 



dx. 



X = R(e, Z)Z 1 / 3 b/a is the solution of Xo(x)/x — 
~ea a /bZ 4 / 3 e 2 , hence X = X{ea a /Z^ 3 e 2 ). We find that 
3tf(c, Z) obeys the following scaling rule, 



g TF (e,Z) = Z-^°± f 



ea 



(26) 



where / stands for the integral appearing in (|25j). 

Consider the single particle density of states at en- 
ergies comparable to the ionization energy, e ~ ~Ei oril 
since only electrons in this energy range participate in 
real scattering processes contributing to the e-e lifetime 
of excitations below the ionization threshold. We refer 
to such electrons as active electrons. In order to evalu- 
ate the TF single particle density of states, Qtf{^i Z), for 
e ~ — e 2 /ao, we divide the spatial integration in (24) into 
two parts: 

a. Core (r < 10 Z _1 / 3 ao) - In this region Vtf(^) ~ 
Ze 2 /r, and |e| <C |VtfM|- The contribution of 
this region to gxF can thus be approximated by, 
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, , (2m) 3 / 2 l-UK-'"". ff^ 
7 -i/3 £o 

The physical meaning of this is that for large Z the 
core electrons lie deep in the atomic potential well, 
and their contribution to the single particle density 
of states at high energies is negligible. 

b. Outer shell (10 Z" 1 / 3 a < r) - Substituting the 
asymptotic expression for the potential, Vtf{t) — 
— 144e 2 a 3 l /& 3 r 4 , valid for r 3> Z _1 / 3 a , we obtain, 



, , 144e 2 a 3 l , . . 
r 2 \U+ - , ° dr (28) 



b 3 r 4 



50 



a I e 



2 \ 1/4 



\eao 



Thus, the TF model indicates that most electrons be- 
long to the core, and only a small number of them (<C Z) 
are active. For a neutral atom the number of active elec- 
trons is estimated by, 



2 /2a 



g T F(e,Z)de. 



(29) 



The contribution of the core electrons to NjfAZ) is 



Z 1 / 3 , so for large atoms Nj^ is composed solely of 



electrons from the outer shell. 



Nei'f'(Z) : 




FIG. 1. Dotted line: Nj f F f (Z) as a function of Z (eq||). 
Solid line: 2-Z 1 ^. Triangles: N e ff(Z) as calculated from the 
discrete single particle spectrum. 

The results of numerical calculation of N^FJZ) are 
shown in figure p), for Z values characteristic of real 
atoms. Indeed Nj f F f (Z) < Z. In this range (Z < 100), 
N^K is found to be roughly proportional to Z 1 / 3 . How- 
ever, within the TF model, NjK(Z) does not diverge as 
oo. Rather it is upper bounded (ill by, 



Z 



M Z N^f f (Z)< 40. 



(30) 



For large enough Z the number of active electrons in 
a TF atom is constant, independent of Z. There are, 



however, no real atoms with large enough Z to check 
this prediction experimentally. 

For a particular atom, the actual number of active elec- 
trons, N e ff(Z), depends irregularly on Z due to the shell 
structure. The number of electrons in the ground state 
configuration with single particle energy e > — e 2 /2ao, ac- 
cording to a Hartree-Fock-Slater calculation |l2| , is com- 
pared with N^f f (Z) in figure (g). It is N eff (Z) and 
not Njff(Z), which determines the spectra of different 
atoms. Nevertheless, NjF(Z) gives a rough estimate for 
N e ff(Z), and leads us to the important observation that 
this number is <C Z. 



2. The Spectrum of Ho 

The TF model suggests that the single particle den- 
sity of states of an atom (or ion) is finite at the atomic 
Fermi energy, even for Z — > oo. However, because of 
the degeneracy of the single particle spectrum and the 
importance of e-e interaction, the characteristics of the 
many-body spectrum vary strongly between atoms. We 
recall that the density of states of Hq is a convolution 
of the density of states of the excited electron, g e , and 
the ionic many-body density of states, gi on - Apart from 
the few lowest levels, the spectrum of the excited elec- 
tron is hydrogen-like in all atoms. Hence, g e depends 
only weakly on the specific atom under consideration. 
The converse is true for the ionic density of states since 
the many-body density of states in the vicinity of the 
ground state strongly depends on the number of valence 
electrons and on the number of available single particle 
states of similar energy. We discuss two limiting cases - 
alkali atoms characterized by a sparse spectrum in the 
vicinity of the ground state, and complex atoms charac- 
terized by an exponentially large density of states in that 
energy range. 

Alkali atoms are characterized by a single valence elec- 
tron in an s-shell. All other electrons reside in closed 
shells with relatively deep single particle energies. Ex- 
citing an additional electron has a relatively large energy 
cost, hence the spectra of singly ionized ions have large 
gaps (nearly ~ Ei) at the ground state. 

Atoms with a few valence electrons are termed "com- 
plex atoms" . Typically such atoms have several open 
shells with similar single particle energies in the vicin- 
ity of the Fermi level. As Z increases these shells can 
acquire large orbital momenta, and consequently large 
degeneracy. The distribution of N e ff — 1 valence elec- 
trons of the ion among those iV s single particle states 
gives rise to ~ N s c/ ' many-body states of comparable 
energy, leading to an average many-body level spacing, 

A ~ jV7 (JVe//-1) e 2 /ao- 

The effect of e-e interaction on the many-body spec- 
trum of a complex atom is demonstrated in the numerical 
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calculation of the spectrum and eigenstates of a Cerium 
atom by Flambaum et al. A Cerium atom contains 
only 4 valence electrons, which are sufficient to gener- 
ate a complex many-body spectrum. The eigenstates of 
a Cerium atom are shown to become chaotic superposi- 
tions of Slater determinants belonging to different config- 
urations, already at low excitation energies. It is evident 
that e-e interaction can not be treated perturbatively in 
complex atoms. Generally (a) it lifts most degeneracies 
peculiar to the non-interacting spectrum, leading to a 
more homogeneous spectrum and (b) it mixes Slater de- 
terminants, so that exact eigenstates are superpositions 
of many determinants. The single particle selection rules 
are hence relaxed, and only rules concerning the total 
spin and angular momentum hold. Consequently, the 
density of states for transitions is significantly increased. 

B. The Interaction Matrix Elements 

We consider matrix elements of Hint between the ini- 
tial excited state, \%) = \n,i, U, rrii) \cti), and any final state, 
|/) = \nf,lf,mf}\af). \ni,k,mi) and |n/,Z/,m/) refer, 
respectively, to the initial (hydrogen-like) and final state 
of the excited electron p3[ . Similarly, |otj) and \otf) refer 
to the initial and final state of the ion. \oti) and \oif) 
are exact eigenstates of the ion, so that \i) and |/) are 
eigenstates of H . 

The multipole expansion of Hint is, 

^ = E E l drl w+i 6 ^'^ (3D 

k— mk— — k > 

xY k , mk (6, 4>)Yl mk {e>^'). 

where Yfc, mfc are the spherical harmonics, r < = mm(r, r 1 ) 
and r> = max(r,r'). Since the excited electron is prac- 
tically exterior to the ion, r> can be identified with its 
coordinate and r< with the ion electrons. The interaction 
Hamiltonian is hence a sum of terms acting separately on 
the excited electron and on the ion, 

A:— rrik — — k 

x J d?r' k 8 P {?)Y^ mh {e',4>'). 

Eq. ( |32| ) is essentially an expansion in the ratio be- 
tween the characteristic radius of Sp(r) (~ do) and the av- 
erage radius of the excited electron (~ n 2 ao). This ratio 
is small due to the large spatial extent of the excited elec- 
tron, which follows directly from the softness of the con- 
fining atomic potential at large distances. Roughly, we 
expect the k th term in the sum to scale like ~ l/n 2t ^ k+1 \ 
Thus, high order terms in the multipole expansion dimin- 
ish quickly. 



Calculation of matrix elements of the simpli- 
fied form of Hint (^) involves separate calcu- 
lations of matrix elements for the excited elec- 
tron, (ni,/i,TOi|p^rrYfe, mfc |n./,Z / ,rn/), and for the ion, 

/ dr'r' k Y£ mh (9 , ,<j>'){a i \6p(r')\af). The infinite sums in 
( |32| ) practically contain only few terms, due to angular 
momenta selection rules. The sum over mk reduces to a 
single term with mk = nif — rrii. As a result of parity 
conservation, the sum over k has non- vanishing terms for 
either odd or even values of k (but not for both). The 
k = term (monopole) vanishes since it is already in- 
cluded in Ho. Due to angular momenta addition rules 
for the excited electron, k is upper bounded by U +lf. 
The terms in the expansion decrease with k, so the dipole 
term is dominant (unless it vanishes). 

Next, we examine the extent to which the interaction 
couples degenerate eigenstates of Ho. The odd terms in 
the multipole expansion, and in particular the dominant 
dipole term (k = 1), do not couple degenerate states. 
This follows from the observation that terms with odd 
k contribute only when the parity of \oii) and \aj) are 
different, and ionic states of opposite parity are non- 
degenerate (apart from rare cases of accidental degen- 
eracy). Hi n t may have non- vanishing matrix elements 
within a degenerate sub-space originating from the even 
k terms in the multipole expansion, prominently the 
quadrupole term (k = 2). In these cases, one should 
use degenerate perturbation theory. However, since the 
quadrupole matrix elements are negligible compared with 
the dipole contribution, we focus on the latter using non- 
degenerate perturbation theory. 

Consider the matrix elements of the excited electron. 
In cases where the final state can be described as an 
ion and a hydrogenic electron, just as we assume for the 
initial state, the matrix clement is calculated by sub- 
stituting the appropriate hydrogenic wavefunctions and 
performing the integrals. This is done in appendix |A[ As 
expected, the matrix elements decrease with increasing 
k, roughly as ~ n~ 2 ( k+1 \ Moreover, the matrix elements 
decrease rapidly as the difference between n, and n f in- 
creases, because the radial integral is strongly suppressed 
by the phase difference between the initial and final hy- 
drogenic wavefunctions. 

When the final state, |n/,Z/,m/), is not a hydrogen- 
like state, the calculation of the matrix element depends 
on the specific atom and particularly on the potential 
at r ~ ao- However, for such final states, the matrix 
elements are negligible because of the small spatial over- 
lap between the confined final wavefunction and the rel- 
atively extended initial wavefunction |l4| ]. 

To examine the dominant dipole term (k = 1) we ex- 
ploit the commutation relation, [h e , Yi, m ]- On one hand, 

(ni,k, Tni\[h e ,Yi tm ]\nf,lf,mf) (33) 
= i e i - e f) (ni,k,mi\Yi tm \nf,lf,mf), 
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while on the other hand, 



(ni,k, mi\[h e ,Yi >m ]\nf,lf,mf) 



(34) 



1*2 



(ni,k,mi\[-^,Y ltm }\nf,lf,mf) 



2m 

2 k(li + l)-l f (l f + l) 

e °o 7T^ 



x (m, ^,mj|-^|n/, //, to/), 



leading to, 
{ni,k, mj|— ^|n/, //, to/) 



(35) 



_ 2 ej-e/ 

e 2 a 



Thus, the dipole term vanishes unless If — k ±1. For 
hydrogen-like wavcfunctions we find, 



(riiJi, TOj|— lj p-\n f , k±l,m- rrn) 



1 1 

a 2 n% 



(k,mi\Yi tm \li ±l,m— rrii) 



(36) 
(37) 



x J dru nztU {r)u nft i i± i(r). 

The overlap integral and hence the matrix element is 
strongly suppressed when n/ differs from rij. 

An exact calculation of the ionic part of the interac- 
tion matrix elements is notably more difficult, since it 
depends on the ionic many-body wavefunctions. Due to 
parity conservation, the dipole term in ( f32] ) has only non- 
diagonal matrix elements, 



dfrY^ im {e,4>)(ai\Sp(f)\a f ) (38) 
= (a t \JdrrY 1 * m (6,<j))p(r)\a f ) = (ai\d l m \a } ). 

where {<i„j}m=±i,o is the tensorial representation of the 
dipole operator, d = J drrp(r). From the Wigner-Eckart 
theorem we learn that d connects only states with equal S 
and M s , AX = ±1, 0, AM L = ±1, and opposite parity. 

We now estimate the dipole matrix elements 
|(ai|d|a/)| 2 for final ionic states that appear in the golden 
rule sum which is restricted by energy conservation. The 
energy transferred by the excited electron ( < Ei) is too 
small for exciting core electrons. Therefore, the final ionic 
states that should be summed over involve valence elec- 
trons excitations only. 

The dipole operator can be expressed in terms of single 
particle creation and annihilation operators, 



dij c \ c j, 



(39) 



where dij are the matrix elements of the single particle 
dipole operator. This sum can be divided as follows, 



E 4 



C Ci 



E d 



ijCiCr (40) 



i,jEopen shells 



i or j£core 



We define a projected dipole operator, d pro j — 
dijc\cj. For final ionic states with no core 

i,j£open shells 

excitations, |a/), 



(a f \d\ 



a i) = \&f\ a proj\&i), 



and, 



E K a /l<W/K)l 2 = («il4r- OJ -K 



(41) 



(42) 



The last equation can be used to estimate the average 
squared dipole matrix element appearing in the golden 
rule, 



|(a/|d]e 



((*i\d%roj\®i 

N tnf . 



(43) 



where Ntot, as defined earlier, is the number of final ionic 
states coupled by the interaction to the initial state. 

Similarly, we define a projected density operator, 
Pproj = c l c j- With this definition, 

i,j£open shells 



d proj = J dr dr' r ■ r 1 p vro j{r)p prol {r') (44) 
df f 2 Pp ro j (r) + / dr dr' f ■ r' g(r,r'), 



(ai\d proj \ai) = J dfr (aj|p proi (r)|a!i) 

+ / drdr'r-r'(ai\g(r,r')\a.i). 



where we have introduced the pair correlation operator 
g{r,r') = p pr0 j(r) Pproj, (r') - 8{r - r')p proj (r). 
Thus. 



(45) 



The first integral in (|46|) is due to auto-correlation while 
the second integral reflects correlations between elec- 
trons. For a system of many electrons interacting re- 
pulsively, one expects the ground state to be a Fermi- 
liquid-like state, with a short range pair correlation func- 
tion g(r, f") (on the order of the average distance be- 
tween electrons). Furthermore, since repulsive interac- 
tion tends to keep the valence electrons away from one 
another, g(f,r^) < for small \r — f |. Consequently, 

{ai\d^ roj \ai) < (ai\ J drr 2 p proj (f)\ai). (46) 

Alternatively, (ai|(2 roj -|a,) can be estimated using the 
chaotic nature of the many-body eigenstates. \cti) can 
be decomposed into a sum of single Slater determinants, 
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\ a i) = Em A ml<W> where A m = (S m \ai). Following 
Flambaum [|| we assume that \cti) are chaotic superpo- 
sitions, so the coefficients A m are random. The number 
of principle components in the decomposition is denoted 
by TV, so that |A„ 



■7JT- Thus ' 



{ai\d proj \oii 



= Y,A* m A n {5 m \d? proj \5 n ) (47) 

m,n 

= ^\A m \ 2 (6 m \d? proj \6 m ) 

m 

It is easily shown that the matrix elements of dp^ ro j in a 
basis of single Slater determinants are all 0(N e ff), where 
N e fjis the number of valence electrons. The first sum 
in (pTj) is 0(N e ff), since it has Af non-negative terms of 
order jjO{N e f /). The second sum is of the same order, 
since it is a random sum of M 2 terms of order jjO(N e ff). 
Both points of view lead, then, to the estimate, 



\(<Xf\d\<XiW 



NeffR 



vol 



(48) 



where R va i is the average radius of a valence electron. 
The number of valence electrons, as well as their average 
radius, vary between atoms due to the shell structure. 

A very rough estimate of the interaction matrix ele- 
ment, taking into account only the dominant dipole con- 
tribution (for |i) and |/) satisfying all the selection rules 
described above), is obtained by combining ( |36| ) and (48), 



\mint\f)f 



1 



Nt. 



2 



e 4 N eff R 



(49) 



N eff R val 

a J 10(nfk)2 : 



for rif close to rij. The dependence upon the specific 
atom under consideration comes from the ionic part of 
the matrix element via N to t , N e f / and R va i . The de- 
pendence of the matrix elements on the initial excitation 
energy is mainly due to the excited electron. Since N to t 
grows exponentially with N e ff, the interaction matrix el- 
ements are very small compared with the level spacing of 
the single particle Hamiltonian, h e . 



excitations in this energy range is very large (infinite in 
principle), because the single particle spectrum becomes 
hydrogen-like as the ionization threshold is approached. 
The existence of many single particle excitations with 
no e-e broadening in alkali atoms follows from the non- 
uniform single particle density of states, characterized by 
a gap at the Fermi energy and becoming dense near the 
ionization threshold. 





1st 2nd 3rd 






generation generation generation * * 


ionization 
threshold 


pfj 






u 






a 

u 


g.s. 





FIG. 2. A schematic illustration of the hierarchical tree for 
an alkali atom. 

As mentioned earlier, the problem of e-e lifetime of sin- 
gle particle excitations can be mapped onto the problem 
of localization on a hierarchical tree. Let us consider the 
hierarchal tree of an alkali atom shown schematically in 
figure (Q). Any n-particle n-hole excitation with n > 1 
involves excitation of n — 1 electrons from closed shells, 
and therefore a large excitation energy. This is mani- 
fested in the tree by a large gap between the lowest state 
in the n th generation and the ground state, a gap that 
grows rapidly with n (see fig. ||). The structure of the 
hierarchical tree forces low lying single particle excita- 
tion to be localized in Fock space. Therefore, they are 
not broadened by e-e interaction. The suppression of 
e-e broadening is characteristic of systems with a gap 
at the Fermi level, and it persists to excitation energies 
considerably larger than the gap. The atomic system 
is unique in that the single particle spectrum becomes 
dense above the energy gap, so the number of single par- 
ticle excitations in alkali atoms for which e-e broadening 
is suppressed is very large. 



2. The Complex Atom 



C. e-e Lifetime 

1. Alkali Atom 

The considerations above show that in alkali atoms 
below the ionization threshold there is no level broad- 
ening due to e-e interaction, and each low lying many- 
body eigenstate can be identified with a particular single- 
particle excitation (i.e. there is a large overlap *~ 1 be- 
tween these two states). The number of single particle 



We argued earlier that the density of states of ionic 
excitations in complex atoms is exponentially large. The 
question is whether this density of states is sufficient 
to induce broadening of single particle excitations. The 
general criteria for applying the golden rule to the cal- 
culation of lifetime is that the resulting scattering rate, 

2 

h/i~ee ^ | (-Hint) | G final (E) (where I (H int ) I is the average 
interaction matrix element and G final {E) is the relevant 
density of states) is larger than the level spacing of the 
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relevant final unperturbed states 1/G finai(E). This re- 
quirement is equivalent to the condition, 



1 



Gfinal(E) 



< \(H int )\. 



(50) 



Gfi n ai(E) refers only to final states which are coupled 
to the initial state by the interaction, and generally may 
be considerably smaller than the full many-body den- 
sity of states. In our case, Gfi na i{E) is smaller than 
the full many-body density of states due to the selection 
rules associated with symmetry (for the dominant dipole 
term). Nevertheless, Gf ina i(E) still grows exponentially 
with the number of active electrons N e ft. N to t is the 
number of accessible final states, and the relevant den- 
sity of states is roughly given by G fi na i(E) ~ N tot a /e 2 . 
We saw in the previous section that \(H int ) \ cx l/^/N to t- 
Since 1/Gfj na i (E) cx l/N to t, eventually for large enough 
Ntot, ( |50| ) is satisfied. We conclude that for large enough 
Ntot, e-e interaction leads to finite lifetime of single par- 
ticle excitations. 

The lifetime of a single particle excitation \i) with en- 
ergy E is, 



1 _ 2tt ^ 



9ion(E + 



2aonj 



)\mH mt \f)\ 2 - (5i) 



The summation is restricted to final states of the excited 
electron, as we have averaged the interaction matrix el- 
ements over the ionic states and correspondingly intro- 
duced the ionic density of states. Taking into account 
only dipole matrix elements (eq. 15^), one obtains, 



1 



27r / N tot ao 

n v e 2 

2n e 2 N eff k Za 
h a lO(nfZi) 5 



e 4 N eff R. 



2 

vol 



(52) 



The single particle excitation spectrum is hydrogen-like 



with E„ 



-/2aonf, so the single particle level spac- 



ing is cx l/n 3 . For large enough n^, the level spacing of 
the excited electron becomes smaller than the level spac- 
ing of the ion. For such single particle excitations, no e-e 
lifetime can be defined and ( |5l|) becomes irrelevant. 

The calculated e-e scattering rate depends on the spe- 
cific atom through N e ffR^ al , which depends irregularly 
on Z due to the shell structure. For all atoms (Z < 105), 
the number of valence electrons < 14, and the resulting e- 
e scattering rate is much smaller than the single particle 
level spacing. 

One may consider atoms with larger Z . According to 
the TF model (see |lIIAl| above), N eff cx Z 1 ' 3 for small 
Z, but in the limit Z — > oo it saturates to a constant. 
This yields e-e scattering rates which are independent of 
Z in the limit Z — > oo. The TF model, however, ignores 
the shell structure. The shell structure suggests that the 



largest open sub-shell in an atom may have Z 1 / 3 elec- 
trons, so the number of valence electrons fluctuates, as 
a function of Z, between 0(1) and 0(Z 1 ^ 3 ) electrons. 
The average radius of the valence electrons, R va i, also 
depends on Z. As the number of valence electrons in- 
creases, we expect R va i to decrease because the screening 
of the nucleus is less efficient. This tends to reduce the 
fluctuations in N e ffR^ al as a function of Z. 

Thus, e-e scattering rates in atoms fluctuate with Z. A 
small fraction of atoms have 0(Z 1 ^ 3 ) valence electrons. 
For such atoms the typical scattering rates increase as a 
function of Z (l/r ee cx Z a , a < 1/3), and eventually, for 
large enough Z, become larger than the single particle 
level spacing. For most atoms, however, the number of 
valence electrons is small and the e-e scattering rates of 
single particle excitations remain smaller than the single 
particle level spacing. 

For a specific complex atom, the width of single par- 
ticle excitations decreases as the excitation energy in- 
creases, in contrast with the analogous width in quan- 
tum dots. This decrease is a direct consequence of the 
soft atomic potential at large distances (r>oo). A small 
increase in the excitation energy amounts to a substan- 
tial increase of the average radius of the excited electron. 
As a result, the interaction matrix elements and therefore 
the e-e broadening are significantly reduced. 

IV. CONCLUSIONS 

We estimated the e-e scattering rate of single particle 
excitations in atoms below the ionization threshold. We 
conclude that e-e interaction does not limit the number of 
discrete single particle excitations observable below the 
ionization threshold in naturally occurring atoms, since 
the broadening it induces is small compared with the 
level spacing of single particle excitations. In practice, 
broadening due to electron-photon interaction limits the 
number of resolved single particle excitations. We find 
that e-e scattering rate in atoms decreases with excita- 
tion energy (below the ionization threshold), contrary to 
its behavior in Fermi liquids. 

The e-e scattering rates in ballistic (clean) and diffu- 
sive (dirty) QDs have been calculated theoretically M. 
These calculations indicate that e-e interaction severely 
limits the number of observable levels in QDs. In a bal- 
listic QD containing N electrons, only ~ y/N discrete 
single particle resonances can be resolved. The effect of 
e-e interaction is even more dramatic in diffusive QDs. 
The number of resolved single particle resonances in a 
0-D diffusive QD Eq] is approximately equal to the di- 
mensionless conductance g which is typically *~ 10 — 20. 
In addition, the e-e scattering rate of single particle ex- 
citations in QDs is shown to increase as the excitation 
energy becomes larger. This Fermi-liquid like behavior is 
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in contrast with our results for atoms. The angular integral is expressed in terms of Wigner's 3-j 

The profound difference of the e-e scattering rate in symbol, 
QDs and in atoms can be attributed to the different con- „ 

fining potentials. In atoms, the confining potential is very / duY l * rn ,(0,(j))Y kimk (0,(j))Y^ m (9,<p) (A2) 

soft at large distances ~ 1/r, while in QDs, the potential J 

has sharper boundaries ~ r 2 . The characteristics of the / {21' + l)(2fe + 1) (2/ + 1) 

confining potential determine the spatial extent of the V 47r 

single particle wavefunctions as well as the single parti- f V k I \ f V k I 

cle spectrum. These properties have substantial influence x I o / \ —to' m k m 

on e-e scattering rates. 

In QDs the confining potential generates a fairly uni- The radial integral gives, 
form single particle spectrum in the vicinity of the Fermi 

/*°° 1 

level. The single particle spectrum leads to a relevant / dr k+1 u n ',i'(r)u n ,i(r) (A3) 

density of states (in this case, the 2-electron 1-hole den- 

sity of states) which grows quadratically with excitation f°° 1 / (n — I — 1)! _ r /„ ao / 2r \ +1 

energy. The resulting e-e scattering rates increase with — J Q r k+i v aQn 2( n _|_ ni \ na o/ 

excitation energy. The spectrum consists of a small num- , 

ber of discrete single particle resonances, beyond which x ^2i+i (2 r /na ) / — ~ ^ ~ ll: c -r/n'a 

the broadening exceeds the single particle level spacing, n-l-i y aQn '2^ n i _|_ 

and the single particle resonances merge and form a con- { r +1 
tinuous spectrum. x ( ) L^'+l^r/n'ao) 

In contrast, the Coulomb confining potential in atoms 
produces a dense single particle excitation spectrum as 
the ionization threshold is approached. The spatial = 
extent of an excited electron's wavefunction increases 
rapidly with excitation energy, due to the softness of the 
Coulomb potential at large distances. The interaction 
matrix elements hence decrease with excitation energy 
leading to reduced e-e scattering rates. The number of 
observed discrete levels is found not to be limited by e-e 
interaction. 
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APPENDIX A: 

The matrix element of Yk,m k (0, (f>)/r k+1 between two 
hydrogenic wavefunctions |n, I, m) and \n' , V , to') encoun- 
tered in section ( |IIIB| ) is calculated in the following way. 
The radial and angular integrations are separated, 

(n', I', m'\Y k>mk (9, 0)-^|n, /, to) (Al) 

,oo 1 

= J dr ^T u n',i'{r)u n ,i{r) 

x f dLuY^ m ,(9,4>)Y kimk (e,0)Y lim (0^). 
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